home *** CD-ROM | disk | FTP | other *** search
/ Linux Cubed Series 3: Developer Tools / Linux Cubed Series 3 - Developer Tools.iso / devel / lang / lisp / gcl-1.000 / gcl-1 / gcl-1.0 / mp / mp_mulul3.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-05-07  |  1.7 KB  |  70 lines
  1.  
  2. /*          Copyright (C) 1994 W. Schelter
  3.  
  4. This file is part of GNU Common Lisp, herein referred to as GCL
  5.  
  6. GCL is free software; you can redistribute it and/or modify it under
  7. the terms of the GNU LIBRARY GENERAL PUBLIC LICENSE as published by
  8. the Free Software Foundation; either version 2, or (at your option)
  9. any later version.
  10.  
  11. GCL is distributed in the hope that it will be useful, but WITHOUT
  12. ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  13. FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  14. for more details.
  15.  
  16. You should have received a copy of the GNU library general public
  17. license along with GCL; see the file COPYING.  If not, write to the
  18. Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  19. */
  20.  
  21.  
  22.  
  23. #include "include.h"
  24. #include "arith.h"
  25.  
  26. /* ulong a,b,y;
  27.    (y = mulul3(a,b,&hiremainder), hiremainder:y == a*b) is TRUE.
  28. */   
  29. #ifdef USE_WORD_MULUL3
  30.  
  31. int mulul3(x,y,hiremainder)
  32.      ulong x,y,*hiremainder;
  33. {
  34.   ulong xlo,xhi,ylo,yhi;
  35.   ulong z; TEMPVARS
  36.  
  37.   xlo=x&65535;xhi=x>>16;ylo=y&65535;yhi=y>>16;
  38.   z=addll(xlo*yhi,xhi*ylo);
  39.   *hiremainder=(overflow)?xhi*yhi+65536+(z>>16):xhi*yhi+(z>>16);
  40.   z=addll(xlo*ylo,(z<<16));*hiremainder+=overflow;
  41.   return z;
  42. }
  43.  
  44. #else
  45. ulong
  46. mulul3(a,b,h)
  47. unsigned int a,b, *h;
  48. {unsigned int temph,templ,ah,al,i;
  49.  ah=0;
  50.  al=0;
  51.  /* in case the shift by 32 does not zero an unsigned int..
  52.     we separate out the first step.*/
  53.  {if (b & 1)
  54.     {temph=0;templ=a;
  55.      lladd(temph,templ,ah,al);}
  56.     /*    printf("\n%d b=%d a=%d (%d:%d)",i,b,a,ah,al); */
  57.     b=b>>1;
  58.   }
  59.  i=1;
  60.  while(b)
  61.    {if (b & 1)
  62.       {llshift(a,i,temph,templ);
  63.        lladd(temph,templ,ah,al);}
  64.       i++;b=b>>1;
  65.     }
  66.  *h=ah;
  67.  return al;
  68. }
  69. #endif
  70.